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Abstract 

We investigate the bounds on the mass of the lightest neutralino from rare meson decays within 
the MSSM with and without minimal flavor violation. We present exphcit formulae for the two-body 
decays of mesons into light neutralinos and perform the first complete calculation of the loop-induced 
decays of kaons to pions and light neutralinos and B mesons to kaons and light neutralinos. We find 
that the supersymmetric branching ratios are strongly suppressed within the MSSM with minimal 
flavor violation, and that no bounds on the neutralino mass can be inferred from experimental data, 
i.e. a massless neutralino is allowed. The branching ratios for kaon and B meson decays into light 
neutralinos may, however, be enhanced when one allows for non-minimal flavor violation. We find new 
constraints on the MSSM parameter space for such scenarios and discuss prospects for future kaon and 
B meson experiments. Finally, we comment on the search for light neutralinos in monojet signatures 
at the Tevatron and at the LHC. 
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I. INTRODUCTION 

Supersymmetry (SUSY) is an attractive candidate for physics beyond the standard model 
(SM) [1-3]. In the minimal supersymmetric extension of the SM (MSSM) the superpartners of 
the electroweak gauge and Higgs bosons (the gauginos and higgsinos) mix to form two electrically 
charged and four electrically neutral mass eigenstates, called charginos, Xi^2' neutralinos, 
Xi 2 3 45 respectively. In many supersymmetric models, the neutralino Xi is the lightest super- 
symmetric particle (LSP) and thus plays a special role in phenomenology. In particular, in 
models with conserved proton hexality [4] (or conserved i?-parity [5]) a neutralino LSP is stable 
and provides a promising dark matter candidate [6-8]. 

So far, no evidence for supersymmetric particles has been found, and lower limits on sparticle 
masses have been derived from searches at colliders and from precision measurements of low- 
energy observables. The particle data group (PDG) [9] quotes a lower limit on the mass of the 
lightest neutralino [10] 

m^o > 46 GeV (95% C.L.). (1) 

This limit, however, is not model independent but assumes the unification of the gaugino mass 
parameters Mi and M2 at some high energy scale. The renormalization group evolution of the 
mass parameters then implies 

Ml = l^M2 = I tan^ ^^Ms ^ ^Ms (2) 

at the electroweak scale, where g and g' denote the SU{2)l and U{1)y gauge couplings, respec- 
tively, and 9w is the weak mixing angle. The experimental mass bound for the lighter chargino, 
m-± > 94 GeV (95% C.L.) [9] places a lower bound on M2 (and on the higgsino mass parameter 
n) and indirectly, through Eq. (2), on Mi. The bounds on Mi, M2 and /i, as well as the lower 
bound on tan/9 = V2/V1 (the ratio of the two vacuum expectation values of the MSSM Higgs 
doublets) from the LEP Higgs searches [11], then in turn give rise to the lower bound on the 
mass of the lightest neutralino, Eq. (1) above. 

In this paper, we investigate a more general MSSM scenario where the assumption Eq. (2) is 
dropped and the gaugino mass parameters Mi and M2 are treated as independent parameters. 
It has been shown that in such a scenario, one can adjust Mi and M2 such that the lightest 
neutralino is massless [12-14]. Specifically, m^o^ = if 

^ _ M2m| sin(2/3) sin^ 9w _ m| sin(2/j) sin^ Ow _ 2 5 GeV ( \ / 150 GeV 
^ /iM2 — m| sin(2/5) cos6'vi/ \tan/5/ \ /i 
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where niz is the Z gauge boson mass. The lower bound on the hghtest chargino mass, 
m-± > 94 GeV, leads to lower bounds of |/i|,M2 > 100 GeV [15]. Thus, adjusting Mi ac- 
cording to Eq. (3) so that m^o = implies Mi -C M2,/i. In this region of parameter space, 
the bino component of the lightest neutralino is above 90%, i.e. the lightest neutralino couples 
predominantly to hypercharge [14, 16]. 

Note that choosing a renormalization scheme in which m^o is an input parameter [17, 18] implies 
that a small Xi mass at tree level will remain small also after radiative corrections. However, 
note that a small or zero neutralino mass can not always be obtained in the presence of non-zero 
complex phases for Mi and fi [18]. (One can always choose a convention where M2 is real by 
absorbing its phase into a redefinition of the gaugino fields.) 

The phenomenology of a very light or massless neutralino has been discussed by many previous 
authors, in particular with respect to its astrophysical, see e.g. Refs. [14, 19-21], and cosmolog- 
ical implications, see e.g. [18, 22-26]. Light neutralinos have also been discussed in the context 
of collider searches [15, 27-32], as well as electroweak precision observables [18, 33] and rare me- 
son decays [18, 34-38]. In a previous paper [18], some of the present authors have presented an 
extensive study of these phenomenological constraints. While a very light or massless neutralino 
cannot provide the cold dark matter content of the universe, it is consistent with all existing 
laboratory constraints and astrophysical and cosmological observations. For a comprehensive 
review of the literature on very light neutralinos, we refer to Ref. [18]. 

In the present work we focus on a specific aspect of the phenomenology of light neutralinos, 
namely rare meson decays to light neutralinos. For more than three decades, rare meson decays 
have been considered a sensitive test for the presence of new physics. In this paper, we investigate 
the tree-level decay of pseudoscalar and vector mesons and the loop-induced decay of kaons to 
pions and light neutralinos and B mesons to kaons and light neutralinos within the MSSM. 
A review and an update of the literature has been presented in a recent accompanying paper 
Ref. [18]. We go beyond Ref. [18] and the previous work in several aspects. First, we provide 
details of the calculations and explicit formulae for the two-body decays of mesons M into light 
neutralinos, M that will allow one to estimate the decay width r(M XiXi) ^-Iso in 

generic supersymmetric theories beyond the MSSM. Moreover, we perform the first complete 
calculation of the loop-induced decays /^"^Tr'XiXi ^iid B" -^K~x\x^i in the MSSM with 
minimal flavor violation (MFV) [39]. We show that these decays are strongly suppressed with 
respect to the SM decays K~ -n^vv and B~ K~vv. We also consider non- minimal flavor 
violation in the MSSM and demonstrate that rare kaon and B meson decays do provide new 
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constraints on the supersymmetric parameter space in the case of hght or massless neutrahnos. 
Finally, we comment on the direct search for light neutralinos in monojet signatures at the 
Tevatron and at the LHC 

We note that this MSSM with extremely light neutralinos is not the only model which allows 
for invisible meson decays into beyond-the-SM (BSM) particles. Other models can have similar 
decays, which may be within reach of near-future experiments. The standard model can be 
simply extended to accommodate light gauge singlet scalars which are dark matter candidates 
[40-42], though it requires some fine-tuning. Rare meson decays to the scalars in this model 
have been calculated in Refs. [43, 44]. These decays in more complex, less fine-tuned models 
[45, 46], in which the light scalars couple to new heavy fermions or a new light f/(l)' gauge 
boson, have been calculated in Refs. [47-49]. Models with extra gauged U{1) symmetries also 
allow for light fermionic dark matter candidates [50, 51], and the rare meson decays into these 
fermions have been calculated in Refs. [48, 49]. The next-to-minimal supersymmetric standard 
model (NMSSM) also allows for light neutralinos, and in the NMSSM they remain a dark 
matter candidate [25]; the NMSSM also admits a light pseudoscalar dark matter candidate. 
The rare meson decays into invisible BSM particles in the NMSSM have been calculated in 
Refs. [44, 52] . Another supersymmetric model is that of Ref . [53] , a model with gauge mediation 
of supersymmetry breaking with hidden sectors, and in these hidden sectors reside candidates 
for light dark matter. Decays of mesons into these dark matter candidates have been calculated 
in Ref. [54]. 

The paper is structured as follows. In Sect. II we present the calculation of the decay of 
pseudoscalar and vector mesons to light neutralinos and compare the results with existing ex- 
perimental bounds. Since the decay amplitudes involve virtual squarks which, in the MSSM, 
are generically heavy, rriq > 100 GeV, the branching ratios to neutralinos are small, well below 
the experimental upper bounds. We provide explicit formulae which allow one to study the 
dependence of the decay rate on the generic supersymmetric masses and mixings. In Sect. Ill 
we discuss the complete calculation of the loop-induced decays K'—^n'x^Xi B~^K~XiXi 
in the MSSM with minimal flavor violation. We present numerical results for various super- 
symmetric scenarios and find that the branching ratios are several orders of magnitude smaller 
than the SM processes with neutrinos instead of neutralinos in the final state. The branching 
ratios for the K^^n^XiXi ^^'^ B^^K~XiXi decays may, however, be greatly enhanced when 
one allows for non-minimal flavor violation in the MSSM, as discussed in Sect. IV. We find 
new constraints on the MSSM parameter space and argue that future experiments in the kaon 
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FIG. 1: Feynmaii diagrams for the decay P/V^XiXi^ omitting diagrams with the Xi^ crossed. 

and B meson sector may be able to probe light neutralinos with masses up to approximately 
1 GeV through rare decays. We comment on direct searches for light neutralinos with a monojet 
signature at hadron colliders in Sect. V and conclude in Sect. VI. The appendix provides details 
on the form factors needed for the calculation of the kaon and B meson decays. 

II. MESON DECAYS INTO LIGHT NEUTRALINOS 

In this section, we consider the decays of the electrically neutral, flavorless mesons M = 
tt", 7], rj', p°, uj, (j), J / ip, T into a pair of light bino-like neutralinos, M— ^-x^Xi- We are thus consid- 
ering the decay of an invisible particle into other invisible particles, analogous to meson decays 
into neutrinos in the SM. Such decays have been studied and searched for experimentally, mak- 
ing clever use of kinematics in decay chains. For example, an invisibly-decaying 7r° can be 
reconstructed from a K'^-^tt^ti^ decay [55], while an invisibly-decaying J/ip may be observed 
by tagging on the invariant mass of the di-pion system in the decay ip{2S)^J/tlj + tttt [56]. 
The decays of pseudoscalar and vector mesons to light photino-like neutralinos have first been 
studied in Refs. [3, 34, 35, 38, 57, 58]. More recent calculations addressing B meson and T 
decays into light new matter in the MSSM and beyond can be found in Refs. [47, 59-61]. We 
comment on the comparison with existing results when discussing the individual decays below. 
Note that while decays with an additional photon in the final state such as T^XiXi + 7 ^i^e 
experimentally easier to detect, they suffer from a further suppression by a factor of a. Since 
the branching ratios we obtain for the two-body decays M^XiXi extremely small anyway, 
we do not pursue the radiative decays M— + 7- 

The tree- level contributions to the decays of pseudoscalar (P) or vector {V) mesons, -P/V^XiXi? 
are mediated by uh/r and c/j^.//? t-channel exchange, as shown in Fig. 1. Here un/R/diL/fi denote 
a left- or right-handed up or down squark of flavor i. 
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A. Decays of pseudoscalar mesons 



We first consider the case 7r°— ^-x^Xi? and then generalize to the other pseudoscalar meson decays. 
Since the branching ratio depends on the ratio m^o/mp in exactly the same way for each 
pseudoscalar meson P and differs only by an overall factor, we plot the dependence once, 
normalized to the peak branching ratio, in Fig. 2. Throughout this paper, to avoid unnecessary 
clutter, we remove the superscripts denoting charges for mesons in subscripts, e.g. we write 
for the mass of the neutral pion, rather than m^o, and we ignore the differences in masses 
between charged and uncharged mesons, so is also used for the mass of the charged pion. 
The decay of a pseudoscalar meson into a fermion-antifermion pair requires at least one of the 
fermions to be massive, due to the required helicity flip. In Fig. 2 we see that the branching 
ratio goes to zero for vanishing neutralino mass. Hence pseudoscalar decays offer no bounds for 
models with massless neutralinos. Instead, the peak branching ratio is reached for m^a/nip = 
1/^6 ^ 0.41. 




m~o/mp 



FIG. 2: The branching ratio for a pseudoscalar meson decay P^XiXi against the ratio m-o/mp 
normalized to the peak branching ratio. 



Likewise, the SM background of P—^UiUi, i = e^n^r^ is extremely small. Its branching ratio is 
proportional to the neutrino mass squared divided by the pion mass squared. Even for the pion 
and relatively heavy neutrinos with my. = 1 eV, this is of the order of 10^^*^. 



1. vr0^x?X? 

We treat the neutral pion as a bound state of valence quarks, given in terms of quark field 
bilinears as {uu — dd)/^/2. Denoting the neutralino momenta by ki and k2, the matrix element 
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for this decay is given by 

M = -{XiXil Xi{i\e\auUlPL + i\e\buU}iPB)u + j^^^^^ 



X 



u{-i\e\auULPR - i\e\buURPL)Xi + d{-i\e\addLPR - \e\ibddRPL)xi k") 



e 



2 



2 
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X (0| f 4- + 4-1 ^7V^ + f 4- + ^ I rfV7'c^k°) - ih^h). (4) 
m^, m.i / \ mi mi ' 



In the last line Uj^o{ki) and v^o{k2) denote the standard 4-component Dirac wave functions (for 
the neutralinos) , whereas u, u, d, d refer to the corresponding quarks. We have approximated the 
squark propagators by ^V^^i^/jj; since the pion mass m^r (and in general all the pseudoscalar 
masses mp) <^mg. In order to obtain the last line we made use of the fact that the pion is a pseu- 
doscalar state, I.e. {0\u-f''u\TT^) = {0\dYd\TT^) = 0. Since the neutralino is a Majorana fermion, 
we also have Xi7^Xi = Xil'^'^Xil + Xil^^Xil = identically. We used Fierz identities [62, 63] 
to transform the products of spinor bilinears to the presented form. 

For purely bino neutralinos, the neutralino-quark-squark couplings are \e\agPL for an outgoing 
ql and \e\bqPR for an outgoing qr, and since the gauge couplings are generation- independent, 
we just write for up-type quarks and ad for down- type quarks, with 

-1 2v/2 -V2 

3y/2cos6w' 3cos6w' 3cos6w 

e is the electromagnetic charge of the electron, such that = Ana. We take a = 1/137 and 
the electroweak mixing angle sin^ 9w = 0.23. Eq. (4) can be seen to agree with Eqs. (2) and 
(3) in Ref. [38], if the neutralino-quark-squark couplings are chosen to be photino-like instead 
of bino-like [i.e. au = bu = 2-\/2/3 and ad = bd = — -\/2/3] and the squark mixing fi^ therein is 
set to zero. 

We use the standard definition of the pseudoscalar decay constant [64, 65] to write 

{0\uri'u\n') = -{0\dn'd\A = -^Vn, (6) 

V 2 

where p-j^ is the momentum of the pion. Eq. (6) defines our convention for the value of the pion 
decay constant F^^ as 131 MeV [9] (as opposed to the other common convention of 91 MeV, 
which incorporates the factor of 1 / -\/2) • 

Summing over final spins, the square of the matrix element is given by 



leading to a partial decay width of 

r.,_,.,„ = i^F>.|.,„,,/i-4f!!«y f 4- + J!^ - 4 - Jf-V . (8) 

If we assume a common squark mass rriq, then we obtain a branching ratio of 

BR(.«^iK) = (4.64xlO-")(^^) (^) ^l-(2.20xl0-)(^). (9) 

using = 135 MeV and F^o = 7.8 eV [9]. 

Since there is destructive interference between the up and down quark contributions, we consider 
an extreme case where the down-type squarks are so heavy that they decouple. This maximizes 
the branching ratio to 

n n ^A /'100GeV\'^/ m^o n2 I / m^o n2 

BR(.«^iW) = (9.30X10-") (^^) (n^) V^-P-20xl0-)(^). (10) 

Setting the up squark masses to 100 GeV, the maximal branching ratio is 1.63xl0~^°, when 
m^o = 55.1 MeV. The current experimental upper bound on the decay of the neutral pion to 
invisible particles is 2.7x10"'^ [55], which is at least three orders of magnitude greater. Hence 
we do not find significant constraints on supersymmetric models with extremely light bino-like 
neutralinos from pion decays. 



2- V^XiXi, rj'^XiXi 

We consider the decays of rj and t]' mesons analogously. We take into account mixing between 
the ?7° and rj^ 5'f/(3)-flavor states, as described by Ref. [66], which leads to: 

{0\uri'u\v) = {Oldr^'dlv) = I (-^cos^sF^s - -Lsin^^F,.^ p^^ ^ ^^I/jX 

(0|s7Vs|r/) = z {^oose^F^s - -^sin^iF,!^ p^^ = iFlp^^ (11) 

and 

(0|n7^^n|r/') = (0|(i7^'ci|r/') = i (^sin^s^^s + i=cos^iF,i^ p^^, ^ tF^^^^^, 

{0\sYjhW) = I (^^sin^gF^s + -Lcos^iF,!^ p';, = iPip';, . (12) 
The phenomenological values of the decay constants and mixing angles are given in Table I. 
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Og 0.37 rad 

9i 0.16 rad 

F^s/F^ 1.28 

F^i/F^ 1.17 

TABLE I: Input parameters for rj and 7]' decay [66]. 
This leads to a partial decay width of 



X 



m (,) 




(13) 



If we again assume a common squark mass m^, then we obtain the branching ratios: 



and 



"9 / 

If we decouple the strange squarks to maximize the branching ratio, we obtain for the rj 

1. /100GeV\^/ TTL^o \2 I / m^o \2 

BRh^x;x;,^(2.63xlO-)(_-_) _ (1.33x10-) (^) . (16) 

We use = 548 MeV, m^/ = 958 MeV, r^,tot = 1-30 keV and r^/,tot = 205 keV, respectively [9]. 
Since in the case of the rj' decay the contributions from each flavor of quark constructively 
interfere, we use the expression with a common squark mass, Eq. (15), to obtain a maximal 
branching ratio. Setting the non-decoupled squark masses to 100 GeV, the maximal branching 
ratios are 7.60x10"^^ for r/^x?X? and 3.83x10"^^ ^'_^^o^o^ when m^o = 223 MeV and 391 
MeV, respectively. The current experimental upper bound on the decay of the 7] to invisible 
particles is 6xl0~^ [9], and for the rj' it is 1.4x10"'^ [9]. Thus for t] and rj' decays we also do 
not flnd signiflcant constraints on models with light bino-like neutralinos. 

B. Decays of vector mesons 

We flrst consider the case of p^^XiXi? and then generalize to the other vector meson decays. 
Since the branching ratio depends on the ratio m^o/mv in exactly the same way for each vector 




0.1 0.2 0.3 0.4 0.5 



m^o/mv 

FIG. 3: The branching ratio for a vector meson decay V—>-XiXi against the ratio m^o/my normahzed 
to the peak branching ratio. 



meson V and differs only by an overall factor, we plot the dependence once, normalized to the 
peak branching ratio, in Fig. 3. We note that there is no helicity suppression in this case, and 
so increasing the x? mass only reduces the available phase space and thus the branching ratio. 
Note that ^Si quarkonium states can only decay into gauginos if parity is violated. As photinos 
and gluinos couple to left- and right-handed squarks with equal strength, the decays ^Si{qq) 
'j'j and ^Si{qq) — gg can only proceed if parity violation is introduced by a difference in the 
left- and right-handed squark masses [57]. Binos, however, couple more strongly to right-handed 
particles than to left-handed particles, i. e. parity is explicitly violated, and quarkonium decays 
into binos are possible also for degenerate left- and right-handed squark masses. 



1- P°^X?X? 

We treat the p° as the vector equivalent of the neutral pion [i.e. the same valence quarks, 
{uu — dd)/V2]. The calculation of the matrix element is similar, except we use 

{0\uYl^u\p^) = {0\dYl^d\p°) = (17) 

due to the CP-properties of the p° meson, and we define 



where is the polarization vector of the p^. The factor Fp we obtain by calculating the SM 
decay width of p^^e~^e~ in this approximation to leading order in a: 

FpO^e+e- = — ^ {QuHu + QdHdf = — , (19) 

6 nip 6 nip 
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where Qu = 2/3, Qd = —1/3 (the quark charges in units of |e|) and Hu = 1/V2, Hd = -1/V2. 
Hu^d are the coefficients of the uu and dd quark field bihnears in the quark model wave function, 
e.g. = {uu — dd) / \/2={HPuu + H^dd), see Ref. [67] for further details. With this generic 
notation, Eq. (19) can easily be adjusted for other vector mesons, e.g. for T, Hb = 1, no other 
H, Qb = -1/3, so Tr^e+e- = {^i^a^ /21){F^/mr) [118]. 
Thus for p°— obtain for the matrix element 



M 



e 
"8 



Wx;(^i)7m7^^x;(^2)(0| 



2 



m. 



U'J^U + 




After squaring, averaging over initial spins, and summing over final spins, we obtain 

2 



96 



( 


«1 ^ 












dR / 





{2krppk2-pp + mlkrk2 - 3m^m|o) 



(20) 



(21) 



leading to the partial decay width 



2 

a 71 



192 ^ ^ 



— r„o ,„+„- 

32 f " 



4 



(3/2) 



1 (3/2) 



m. 



ur 



hi 


al 


) 


2 

m - 

dR 


m - 

dL 








al 




2 ' 

m- 

dR 


2 

dL 



:22) 



If we minimize destructive interference by decoupling the left-handed up squark and the right- 
handed down squark, we obtain 



BR(p°^X?X?) 



fS.OlxlO-^^) 



100 GeVV 



(6.65x10" 



1 Mev) 



(3/2) 



(23) 



usmg nip 



775 MeV, Tpo. 



/Fpo^tot = 4.71x10 ^ [9]. The branching ratio is maximized for 



m^o = 0. Setting the non-decoupled squark masses to 100 GeV, we find a branching ratio of 
8.01x10"^^. We can find no experimental bound at this time for the branching ratio of p°— 
invisible, but note that a branching ratio of 10~^^ is far below the bounds on other branching 
ratios measured for the p° meson [9]. 



2. uj^XiXi, (p^XiXi 

The mixing between the SU (3)-fiavor states amongst the vector mesons is very different from 
the pseudoscalars. Conveniently, it is such that the (p meson is nearly pure ss and the u is 
{uu + dd) I v^, and we treat them as these pure states. 
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This means that the decay uj-^XiXi is analogous to p— taking into account the construc- 
tive rather than destructive interference between the up and down quarks: 



32 



m„ 



21 (3/2) 



m. 



2 



2 

m. 



(24) 



If we take a conservative estimate by decouphng the left-handed up and down squark masses, 
we obtain 

J 1 - (6.53xlO~^W 



BR{u^xlxl) = (1-57 X 10 



100 GeV 



1 MevJ 



(25) 



with = 783 MeV, r^^e+e-/r^,tot = 7.16x10 ^ [9]. This branching ratio is maximized 
for m^o = 0. Setting the non-decoupled squark masses to 100 GeV, we find BR(co'— s>XiXi) = 
1.57x10-1^ 

The (f) case is particularly simple. The partial width is given by 



1] 

32 



(3/2) 



m 



sl 



(26) 



If we minimize destructive interference by decoupling the left-handed strange squark mass, we 
obtain 



BR(0^X?X?) = (7.51 X 10 



-14\ 



100 GeV\' 



m 



SR 



■) 



1 - (3.85x10" 



(3/2) 



(27) 



1 MeV. 

= 2.97x10-^ [9]. Setting 



which is maximized for m^o = 0. Here m,^ = 1.02 GeV, T^p^^+e- /'^(t),tot 
the right-handed strange squark mass to 100 GeV, we find a branching ratio of 7.51 xlO~^^. 
Again, we can find no experimental bound at this time for the branching ratios of (p or uj^ 
invisible, but note that branching ratios of 10"^'^ and 10~^^ are far below the bounds on other 
branching ratios measured for these mesons [9]. 



3. j/i^^mi^^mi 

The squarks associated with the heavier flavors of quarks (charm, bottom and top) are expected 
to mix their left-handed and right-handed components much more than the other flavors [2, 3]. 
We take this into account for the J / tp and T decays. Writing the mass eigenstates of the charm 
squarks Ci^2 as 

'ci\ _ I cosOc sinOc \ f cl 
a) \ sine, cos9, j 



:28i 
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we obtain 



M = ^(XiXil |xiNe|a«(cos6'cC* + smOccDPL + i\e\bu{cos6cC2 - sin6'aCi)PR]c| 
X {c[-i\e\au{cos9cCi + sm9i,C2)PR - z|e|6„(cos6'cC2 - sinOcCi) Pl]Xi} \J/i') 
= —111, 



X 



'Cl 



ml 



(29) 



since 



'\cPl/rc\J/iP) = 0. Hence 



\M\^ = 



48 



g^cosgc)^ - (buSmOcy _^ (a„sin6'c)^ - (&MCos6'g)^ \ ^ 
ml ml J 



x{2krpj/^k2-pj/^ + m]/^krk2 - Sm'^j/^m'io] 



(30) 



and so we obtain 

r 



1-4 



m^o 



J/4> 



(3/2) 



i^cosgg)^ - (5„sin6'g)^ (g^singg)^ - (&^cos6'^ 



2 



(31) 



Our results agrees with the decay width for orthotoponium decay to a photino pair in Ref. [57], 

'u = 2/3, 0-^i^9t, m-ii^mi for the case m^—^Q, and assuming 



with the replacements au/d, bu/d—^e 
mg = 2mt. 

The greatest branching ratio occurs when 6*5 = 7r/2 and the heavier charm squark is decoupled. 
In this limit, and using mj/^ = 3.10 GeV, rj/^^e+e-/rj/i/),tot = 5.94x10"^ [9], the branching 
ratio for J/ip-^XiXi is 



BR(J/V^-.X?X?) = (5. 



(3/2) 



(32) 



U2xlO-)fH^^y[l-(4.17xlO-)(-^ 
\ m-^^ ) \_ Vl MeV 

which is maximized for m-j^o = 0. Setting m^j to 100 GeV, we find a branching ratio of 5.12x 10~^. 
This is well below the experimental upper bound of 5.9x10"^ [69], but not very much below 
the branching ratio for the SM process J/ip-^vv of 2.70x10"^ [61]. 
Analogously, 

(3/2) 



1-4 



mr 



X 



{adCos9i)^ - {bdSinOiY (orfsing^)^ - {bdCos9 
mf mf 

bi 62 



(33) 
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Here 6^ is the mixing angle in the bottom squark sector. The greatest branching ratio occurs 
when 6^ = n j1 and the heavier bottom squark is infinitely heavy. In this limit, the branching 



ratio for T^x'JXi is 



BR(T ^ = (4.47 x lO'^) 



100 GeVV 

1 / 



(4.69x10 



-8n 



1 MeVy 



(3/2) 



(34) 



tot 



which is maximized for mj^o = 0. Here we have used mx = 9.46 GeV, rx^e+e-/rT, 
2.38x10"^ [9]. Setting the lighter bottom squark mass m^^ to 100 GeV, we find a branching 
ratio of 4.47x10^^, which is well below the experimental upper bound of 2.5x10^'^ [70], and 
also well below the branching ratio for the SM process T^vv of 1.05x10^^ [61]. 
We note that T— ^ invisible has been used to probe light dark matter through inversion of the 
annihilation cross-section [47]. 



III. K-^7r-x°x? AND ^-^K-^^X? WITHIN THE MINIMAL FLAVOR VIOLAT- 
ING MSSM 

As mentioned in Sect. II A, the two-body decays of pseudoscalar mesons all have widths propor- 
tional to the square of the mass of the Xi? and so there are no bounds for the case of massless 
neutralinos. This is not the case for three-body final states, and so we examine whether there are 
significant bounds from the decays K^— s>7r~x5Xi aiid 5^— S'-ft'^XiXi! where the SM analogues, 
K^^TT^uu and B^^K'uu, have very low branching ratios [71-74]. 

We present the results for the case of massless neutralinos. The phase space is smaller as the Xi 
mass increases, and any terms in the matrix-element-squared proportional to the Xi mass are 
sub dominant. 



A. General analysis 

For this section, we impose minimal fiavor violation (MFV) [39] at the low energy scale at 
which we are working [i.e. at the mass of the decaying meson). MFV models are defined by the 
requirement that all fiavor violating and CP- violating transitions are described by the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix of the SM, see e.g. Refs. [75, 76] and references therein. 
Thus the decays K~^7r~XiXi ^'^^ B~^K~XiXi only have leading contributions from one-loop 
diagrams. Those for K^-^ir^XiXi are shown in Figs. 4 to 7, without the diagrams with the 
neutralinos crossed. Through the loops, the partial width depends on many parameters, and not 
in a simple way. Therefore we have opted to use a set of supersymmetric benchmark scenarios, 
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thereby fixing the masses and couphngs. We use the masses and mixings of the Snowmass points 
SPSla, SPS2, SPSS, SPS4 and SPSS [77], calculated using SOFTSUSY [78]. We have modified 
the neutralino sector of these benchmark points by choosing Mi according to Eq. (3) such 
that the lightest neutralino is bino-like and massless. The resulting supersymmetric models 
we denote 'pseudo-SPS' points. The imposition of MFV and the degeneracy of the up- and 
charm-squark masses leads to a super-Glashow-Iliopoulos-Maiani (super-GIM) [79] suppression 
in the K~—^7t~XiXi decay and to a lesser extent in the B~^K~XiXi decay. 
Since the scale of the sparticle masses and the W boson mass are so much greater than the 
masses of the kaon or B meson, we integrate out the heavy degrees of freedom to obtain effective 
four-fermion vertices [80]. We then write the amplitude for K~^7r~XiXi 

xixi 2 \ rriK 

+CLHh)PLv{k2){n-\ds\K-) + CRuih)PRvik2){n~\ds\K-)^ - {h^h), 

(35) 
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FIG. 6: Triangle Feynman diagrams for the decay s^dxiX' 



W-/G-/H- 




FIG. 7: Circle Feynman diagrams for the decay s^dxixv 
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recalling that ki and k2 are the momenta of the (indistinguishable) final-state neutralinos. The 
Ci coefficients, i = L, R, are dimensionless. These are the only terms possible, since Xi7^Xi 
and Xi[7^; T'^IXi ^-^e identically zero for Majorana fermions. Conservation of momentum and the 
Dirac equation for the spinors can be used to reduce any other terms to the set above, [applying 
the identities u{k2)'~i^'~i^v{ki) = v'^ {k2)'~i^'~i'^^u^ {k2) = —u{ki)'y'^'y^v{k2) and u{k2)PL/Rv{ki) = 

V^{k2)PL/RU^{k2) = -u{k,)PL,av{k2)]. 

The coefficients C, were calculated with the aid of FormCalc [81], using techniques detailed in 
Ref. [82], with the approximation of neglecting all external momenta compared to the sparticle 
and W boson masses. The loop integrals in this approximation are straightforward and can be 
written as powers and logarithms of masses, so numerical instabilities in the reduction of tensor 
integrals with extreme scale differences are avoided. Variations of these coefficients over the 
phase space of the decay are of order m\- / rn^ and therefore we ignore them. This is analogous 
to the approach taken in Ref. [83] in calculating the leading order value for the standard model 
decay K'^^ir^e'^h'. 

Our calculation takes into account all diagrams at the one-loop level and is not restricted by 
requiring the squark masses to be small compared to the mass of the W boson, thus improving 
the original calculations for photino-like Xi of Refs. [34, 35]. 

The relevant terms for B~^K~XiXi can be obtained from the above by replacing the mesons 
in the bra and ket appropriately and the s quark with a b quark, and the d with a s. 
We use the mesonic form factors as given in Ref. [84] for K'^tt'XiXi R-^f. [85, 86] for 
B~^K~ XiXi (reproduced for convenience in App. A), but for the moment, we shall use the 
shorthand for the mesonic current {'n-^\d'~f^s\K~) = F^. Taking the dot product of both sides 
with the difference of the momenta of the kaon and the pion, which to a good approximation 
is equal to the difference in momenta of the strange quark and the down quark, and using the 
Dirac equation for the quarks allows us to write 

{.-\ds\K-) = (36) 

In this notation 

F-{PK-P7r) 



mKims + rrid) 



As can be seen in Sect. IIIB, the coefficient is much larger than the other coefficients for 
the benchmark points we use for the decay K~ ^Ti~ x\x\- If we neglect all the C coefficients 
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except C^, we find that the spin- averaged matrix-element-squared for K — ^vr XiXi is 

l-M^-V^o^oP = 4G||C^|' iF*-hF-k2 + F*-k2F-h - F*-Fkrh) , (38) 

which differs from that for K'^n^ev (which has a branching ratio of 5.08x 10^^ [9]) only by an 
overall factor [assuming a massless electron and isospin symmetry for the up and down valence 
quarks; the factor of (1/a/2)^ from the coefficient of uu in vr" balances the factor of 1/2 from 
identical final-state particles in the ^tx^ case]. In this case the behavior over the phase 
space is identical and the ratio of the decay widths is equal to the ratio of \Cf\^ to 
For the case of ^K"Xix1, we see that is not as dominant as for K^^n^x^xl, and we 
present the results from directly evaluating the decay width using the expressions for F^ given 
in Ref. [85, 86] (though with diminished accuracy from using mesonic form factors). Direct 
evaluation was found to agree at the percent level for K~—^'7i'~XiXi- For B^^K^xIxi, Cr is 
comparable to and the effect of neglecting it decreases the branching ratio from by roughly 
0.3% for pseudo-SPS5 up to by 66% for pseudo-SPS4. 

B. Numerical results 

We present the numerical values of the C coefficients in Tables II and III. In the final column of 
each table, headed BR/exp, we present the ratio of either the computed K~—^7i~XiXi branching 
ratio to the experimental value of the branching ratio for K~ ^-k^vv (which is 1.73xl0~^° [8^]), 
or the computed branching ratio of B~ -^K~ x\x\ to the current experimental upper bound on 
the branching ratio for B^^K^vv (which is 1.4x10"^ [74]). 



We see that the branching ratios for the decays K~ ^ti' x^Xa. ^^"^ B~ ^K~x\Xa_ iii the MSSM 
with MFV are several orders of magnitude smaller than the respective standard model processes, 
with no chance to observe them in the experiments of today and the foreseeable future. 
The strong relative suppression of the kaon and B meson decays to neutralinos compared to 
their respective SM decays to neutrinos can be explained as a combination of various factors. 
First, the sparticle masses of the SPS benchmark scenarios are considerably larger than the 
W mass; in particular, the SPS2 benchmark point with squark masses of (9(1.5) TeV is most 
strongly suppressed. Moreover, the SM process has four SU (2) vertices while the MSSM process 
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pseudo- 

SPS \Cjf\ \C^\ \Cl\ \Cr\ BR BR/exp 



la 5.70 X 10"^ 7.45 x lO-^"^ 4.93 x 10"^^ 2. 60 x 10"^° 3.28 x 10"^*^ 1.90 x 10"^ 

2 3.81x10-9 7.77x10-1'' 9.55x10-^^ 4.10x10-^^ 1.47x10"!^ 8.49x10"^ 

3 2.63x10-*^ 3.74x10"^'^ 5.13x10"^^ 1.25x10"^° 6.99x10-^'^ 4.04x10-'^ 

4 2.95x10-*^ 6.41x10-1'' 1.34x10-^2 1.42x10-^ 8.76x10-^'^ 5.06x10"'^ 

5 7.12x10-*^ 2.81x10-1"^ 5.61x10-^2 r.UxlO-^^ 5.12x10"^*^ 2.96x10-*^ 

TABLE II: Numerical values for the coefficients Cj defined in Eq. (35) for K~^tt~x\x\ the various 
pseudo-SPS points described in the text. The final column shows the ratio of the branching ratio for 
i^-^vr-XiXi to the experimental value of the branching ratio for K~ ^-k' uu (1.73x10-^^). 

has two SU{2) vertices and two U{1) vertices, resulting in a further suppression of the MSSM 
decay matrix element proportional to tan26'^. Since the flavor-changing aspect is carried by the 
W boson or the charged Higgs bosons, the quarks/squarks are mainly left-handed, and their 
couplings to the bino-like neutralinos are suppressed by their small hypercharge. Finally, the 
spin structure of the MSSM decay amplitude provides another order-of-magnitude suppression 
with respect to the SM decay. Considering all these factors together leads to a suppression of 
several orders of magnitude. 

IV. K- T^~xlxl AND B- K-/tz-xIxI WITHIN THE NON-MINIMAL FLAVOR 
VIOLATING MSSM 

Going beyond the previous section, we here consider non-minimal flavor violation in the MSSM, 
for which there are many potential sources. For example, the squark mass matrix in the soft- 
breaking part of the Lagrangian in general couples squarks of different flavor to each other, 
resulting in flavor changing neutral currents (FCNCs). Thus from the theoretical point of view 
there is a priori no reason to assume MFV within the MSSM. 

To deal with the complicated flavor structure of the MSSM, we employ the mass insertion 
approximation [76, 88], which is defined in the super-CKM basis [76, 88, 89]. In this basis. 
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pseudo- 

SPS \Cj^\ \C^\ \Cl\ \Cr\ BR BR/exp 



la 6.49x10"*^ 3.77x10"*^ 1.93x10"^ 7.22x10^"^ 3.35x10"^° 2.39x10"^ 

2 5.44x10-'^ 3.69x10-9 4.13x10-^ 1.83x10-"^ 2.48x10-^2 

3 3.00x10-^ 1.87x10"^ 2.16x10"'^ 4.54x10-'^ 7.19x10"^^ 5.14x10-^ 

4 3.31x10-^ 3.22x10-9 5.85x10^9 6.38x10-^ 2.53x10"^° 1.81x10-^ 

5 9.50x10-6 1.49x10-9 2.12x10"^ 2.27x10"'^ 7.14x10-^° 5.10x10-^ 

TABLE III: Numerical values for the coefficients Cj defined in Eq. (35) for i^-^il'-XiX? the 
various pseudo-SPS points described in the text. The final column shows the ratio of the branching 
ratio for B~ ^K~^^ to the current experimental upper bound on the branching ratio for ^K~ui> 
(1.4x10-5). 

the Lagrangian is written in the quark mass eigenstate basis. The squarks are rotated with 
the same rotation matrices as the quarks. This leads to a diagonal flavor structure of the 
gauge interactions, but the resulting squarks are not necessarily mass eigenstates. The squark 
propagator can be expanded according to Ref. [90]: 

" (PI - mn - 5m2),^ 

1 % 

= p_^2 ^«/^ + (p_^2)2 ^</3 + OiSm^), (39) 
where 1 is the identity matrix, a and (3 are the indices of the 6x6 squark mass-squared 



matrix and = y {M'^)aa{M^)f3i3 is the "average" squark mass squared. The off-diagonal 
elements of the squared squark mass matrix in the super-CKM basis are given by Sm'^ij. 
To avoid large FCNCs, which would be in contradiction with experimental observations, we 
shall make the common assumption ^ ^"^a/j [76]. In this case higher order terms in Sm"^^ 
in Eq. (39), denoted by 0{6m'^), can usually be neglected. However, if the linear term in 5m^^ 
in Eq. (39) vanishes, higher orders in 6m'^ need to be taken into account. 

In what follows we parametrize the amount of flavor violation with the help of the dimensionless 
mass insertion parameters [76]: 

(5S)xy = ^^%^, ^^J, (40) 
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S Sx {^ds)xY d S 3x i^dJxY d 




FIG. 8: The tree-level diagrams contributing to the process K ir XiXi without assuming MFV. 
The FCNCs are mediated by left- {X, Y = L) and right-handed (X, Y = R) squarks, respectively. 

where (5"^dxy)u with XY = LL, RR, LR, RL are the off-diagonal elements of the 6x6 down- 
squark mass-squared matrices, which couple squarks of helicity X to squarks of helicity Y. 
From the non-observation of FCNCs beyond the standard model contributions, bounds on the 
mass insertion parameters, Eq. (40), can be derived [76, 91-100]: 

\iSi)LL/RR\ < 0(10-') , ms)nL/LR\ < 0(10"=^) , 

\{Si)LL/RR\ < O{10~') , \{6i)nL/LR\ < O{10-') , 

\{Si)LL/RR\ < 0(10-') , \{6i)nL/LR\ < O{10-^) . (41) 
Here an average squark mass of m = 500 GeV is assumed. 

The bounds in Eq. (41) were obtained mainly from meson-antimeson mixing. In this case, 
they scale (at leading order) with m/(500GeV) as long as the ratio of the squark and gluino 
masses is fixed, i.e. the bounds get weaker for increased squark mass. Note that in addition to 
experimental constraints, there exist theoretical constraints on all {Sfj)LR from the requirement 
of vacuum stability [101]. These bounds are comparable or even stronger than those in Eq. (41). 
Furthermore, they do not become weaker as the SUSY scale is raised. 

The constraints on the mass insertions which connect squarks associated with different helicity 
are usually much more restrictive than those on the "helicity-conserving" mass insertions. We 
therefore concentrate in the following on processes where flavor violation is mediated by either 
only left-handed or only right-handed squarks. We thus consider only contributions from either 
Sll or 6rr. 

A. Matrix elements and partial width 

We show in Fig. 8 the tree-level Feynman diagrams which allow for the process K~ t^^XiXi 
assuming a nearly massless bino-like neutralino Xv We obtain, after employing Fierz transfor- 
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mations, the following quark- level matrix elements from the diagrams in Fig. 8: 

= 9cll^etZ ^^"'"^'"'^ |S(«=.)7.T5«(fe)l . (42) 

Here Ai ll and M.rr are the matrix elements involving only left- and only right-handed squarks, 

respectively. The quark currents are replaced by their corresponding hadronic matrix elements, 

which are parametrized by the form factors /+ and /_; see App. A for details. 

From the squared matrix elements, Eq. (A8), we easily obtain the partial width for K~ 

7r~XiXi- Note that we have to introduce an additional factor of 1/2 because this process 

involves two identical particles in the final state. 

The decays 

B- ^-K-lK-xlxl. (43) 

are described by similar squared matrix elements, which can be obtained from Eq. (AS) by 
replacing the form factors according to Eq. (A6) and, if considering B~ -^K~x\}(^, Pn ("^tt) by 
Pk {rriK). 



B. Branching ratios and excluded parameter space 

The SM process K~ 7t~uu has the same experimental signature as the kaon decay into a pion 
plus neutralinos, namely a charged pion and missing energy. Here is a neutrino of arbitrary 
flavor. The theoretical prediction for the branching ratio for this decay (within the SM) is [71] 

BR{K- 7r~i^P)|thcory = (8.5 ± 0.7) x 10"^^ (44) 

The squared matrix elements of the SM decay depend on the external momenta in exactly the 
same way as those of the SUSY decay in Eq. (A8) for massless XiS- The distribution of the 
pion momentum for K~ '^~XiXi is then equal to the p7r-distribution of the SM process. 
However, this is no longer the case for massive x^s as can be seen in Fig. 9. 
We show in Fig. 9 the p,r-distribution [dBK/ dp.j^)/ BK where BR = BK{K~ '^"XiXi)^ fo^^ 
different values of m^o in the kaon rest frame. The distributions are normalized to one. We 
have employed the squared matrix element of Eq. (A7) to calculate ((iBR/(ip^)/BR. We see 
that for m^o > 130 MeV, p-,^ is smaller than 140 MeV. This has important consequences for 
experimental searches. 
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The E787 and E949 collaborations have observed events consistent with the SM decay K — > 
-K^vv. They found that [87] 

BK{K- ^ 7r-i.z>)|exp. = (1.7311:^^) x 10^1°, (45) 

assuming a p,r-spectrum equal to the SM prediction. To separate the signal from the background, 
Ptt regions were selected, namely 211 MeV < < 229 MeV (region I) [73, 102] and 140 MeV 
< < 199 MeV (region II) [87, 103]. These were chosen such that the background K~ tx^tx^ 
with ^ 205 MeV is excluded. It follows from Fig. 9 that the experimental searches in region 
I (region II) were insensitive to the process K~ — > vr"x'j'Xi if "^x? ^ MeV ijn^o > 130 MeV), 
because the respective pion momenta are then too small. 

In the following, we will estimate the experimental sensitivity for scenarios with m^o ^ with 
the help of the correction factor 

(46) 





^0) 




= 0) 




^dBR 


il40MoV 





with 

where Pn,ma.x is the maximal kinematically allowed pion momentum for Xi^ with mass m^o . We 
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FIG. 10: Branching ratios (BRs) for K~ ''^~xixi as a function of the neutralino mass m^o and 
the mass insertion ((^J^^)ll (left-hand figm'e) and ((5^^)_Rij (right-hand figm'e). The branching ratios 
are normahzed to the SM prediction BR(il'~ — > -n^uv) = 8.5 x 10""*^^ [71]. We assume an average 
squark mass of m = 500 GeV. The hghtest gray (yellow) region corresponds to normalized BRs > 10. 
The upper (lower) solid gray (turquoise) line corresponds to the experimentally measured branching 
ratio plus two sigma, Eq. (45), multiplied with the correction factor /c, Eq. (46), (SM prediction) for 
K~ — > 7r~z/i>. The dashed lines show upper bounds on the mass insertions {S'^g)LL/RR obtained from 
mixing [92] for different ratios of squark and gluino masses, cf. Eq. (48). 

set /(m^o) equal to zero if Pn,max < 140 MeV. 

We show in Fig. 10 the branching ratio for K~ 7r~XiXi a function of m^o and the mass 
insertion parameter {6^^)11 (left figure) and {S'^s)rb^ (right figure). We assume in Eq. (A7) an 
average squark mass of m = 500 GeV. The branching ratios are normalized to the SM prediction, 
Eq. (44). 

The lower solid turquoise lines in Fig. 10 correspond to regions in the MSSM parameter space 
where the process K~ — s> 7r~XiXi ^ branching ratio equal to the SM prediction, Eq. (44). The 
upper solid turquoise lines show points in the 'nT'xi-{6'^^) planes, where the branching ratio for 
K~ T^^XiXi is equal to the central experimental value for BR(K~ n^uu) plus 2a, Eq. (45), 
multiplied by the correction factor /c, Eq. (46). We assume, as a conservative approach, that the 
branching ratio for K~ — > ir^uu is negligible compared to the branching ratio for t^^XiXi- 
With our conservative approach, we can exclude at 2a the MSSM parameter space which lies 
above the upper solid turquoise line. As can be seen from Fig. 10, the experiments were unable 
to probe the region with m^o > 130 MeV, although BIi{K~ t^'XaX^) can be one order of 
magnitude larger then the respective SM process, cf. Eq. (44). Future experiments, such as 



24 



CERN P-326/NA62 [104], will be able to explore roughly the MSSM parameter space between 
the two solid turquoise lines as long as m^o < 130 MeV. 

In Ref. [92], bounds on the mass insertion parameters {S'^g)LL/RR were derived from the mass 
splitting of the K^-K^ system. For an average squark mass of m = 500 GeV, three different 
bounds on {S^s)ll/rr were obtained assuming three different ratios x of gluino, rrig, and squark 
masses 



x = ml/m^ = 1.0 =^ (5i)LL//?ii < 4.6 X 10-^ 

x = ml/m' = 0.3 =^ (Si) ll/ rr < 2.2 x 10~^. (45 



The horizontal dashed lines in Fig. 10 correspond to the bounds above. Thus, with our calcula- 
tion of the non-MFV process K~ — ^ t^^XiXi^ we are able for the first time to exclude the MSSM 
parameter space which lies above the upper solid turquoise lines and below one of the dashed 
lines. Note that the weakest bound on {S^g)RR, {x = 4) lies outside the right figure in Fig. 10. 
The plot on the left in Fig. 10 with 7^ looks qualitatively the same as that on the right 

in Fig. 10 with {S'^g)RR 7^ 0. The main difference is that the branching ratio for any specific 
value of (S'^s) and m^o is 16 times larger in the right than in the left figure. This is because the 
light neutralino has to be bino-like and therefore couples to hypercharge. This is larger for right- 
handed particles than the respective left-handed particles in the MSSM. We see in Fig. 10 that 
for a massless Xi? we can strengthen the bound on the flavor- violating mass insertion parameter 
{Sds)RR by a factor of between two and ten depending on the ratio of squark mass to gluino 
mass, cf. Eq. (48). 

As pointed out in Sect. IV A, we can also employ our squared matrix elements, Eq. (A7), to 
calculate the branching ratios of the flavor changing B decays, Eq. (43). Experimentally, so far 
only upper bounds for these processes exist [9, 74, 105]: 

BR{B- ^ n-uu) < 1.0xlO-\ 

BR{B- K^w) < 1.4 X 10^^ . (49) 
The upper bounds on the respective mass insertion parameters are 

x = 4.0 =^ (5i)LL/i?i? < 7.0 X 10-1, (4)LL//?if < 1.1x10°, 

X = 1.0 ^ {5i)LL/RR < 1.4 X 10-\ {5i)LL/RR < 4.8 X 

a: = 0.3 =^ {5i)LL/RR < 2.3 x lO'S 

X = 0.25 =^ {5i)LL/RR < 6.2 X lO-l (50) 
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FIG. 11: Branching ratios for '^~XiXi ^ ^ function of the neutralino mass m^o and the mass 

insertion ((5^f,)LL (left-hand side) and ((5^f,)_R_R (right-hand side). The branching ratio is normalized 
to the calculated BR for B^^tt^uu within the SM, which is 2.2x10^^ [115]. We assume an average 
squark mass of m = 500 GeV. The dashed lines show upper bounds on the mass insertions {5%j)ll/rr 
mainly obtained from B^—B^ mixing for different ratios of squark mass to gluino mass [94], cf. Eq. (50). 

Again, an average squark mass of m = 500 GeV is assumed. The bounds on (5^^) were obtained 
in Ref. [94] from the SUSY contributions to B^-B^ mixing and from the B ^ J/tp Kg CP- 
asymmetry. We calculated the bounds on (5^^) analogously to Ref. [76] from the recently 
measured B^-B^ mass difference [9, 106]. 

Note that additional bounds on (5^^) from Bs Xsi~^i~, Bs — > yLt"'"/i~ and especially from 
b s'-f exist, which can be up to an order of magnitude greater than the bounds in Eq. (50). 
In addition, SUSY contributions to B^-B^ mixing can be dominated by two-loop double Higgs 
penguins for large tan/3. For more details, see for example Refs [76, 95, 99, 100, 107] and 
references therein. However, these bounds are highly model dependent, i.e. they depend not 
only on the squark and gluino masses. We will therefore not consider them in our analysis. 
We show in Fig. 11 the branching ratios for B~ t^'XiXi ^ function of m^o and the mass 
insertions {S^fj)LL (left figure) and {S^i,)rr (right figure). The dashed lines correspond to the 
upper bounds on the the mass insertions as given in Eq. (50). The allowed region for non- 
vanishing {S^f,)LL [{^db)RR]y region below the dashed lines, lies at least two [one] order 
of magnitude below the current experimental upper bound of 1.0 x 10~^, cf. Eq. (49). We can 
thus not exclude additional regions of the MSSM parameter space from the upper bound on the 
branching ratio. However, a future super B factory will be able to explore large parts of the 
parameter space shown in Fig. 11 [108], i.e. branching ratios down to (9(10~^ — 10"'^). Note 
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FIG. 12: Branching ratios for B~ K"x^Xi &s a function of the neutrahno mass m^o and the mass 
insertion {S'^jj)LL (left-hand side) and (5^j)i?_R (right-hand side). The branching ratio is normalized to 
the calculated BR for B~ ^K^uu within the SM, which is 4.5xl0~^ [116, 117]. We assume an average 
squark mass of m = 500 GeV. The solid gray (turquoise) line corresponds to the experimental upper 
bound, Eq. (49), multiplied by a correction factor analogously to Eq. (46). The dashed lines show 
upper bounds on the mass insertions {^'^i^)ll/rr obtained from B^-B^ mixing for different ratios of 
squark mass to gluino mass, c/. Eq. (50). The weakest bound on ((5^^) (x = 4) lies outside the figure. 

that if m^o > 1 GeV, the momentum of the pion or kaon from the B meson decay may be too 
soft to pass experimental cuts. 

Finally we show in Fig. 12 the branching ratios for the process B~ ^ function of 

m^o and (5^^) ll/rr- The numerical values are of the same order of magnitude as those in Fig. 11. 
Differences are mainly due to different form factors, see Eq. (A6). In contrast to the B decay 
into a pion, the decay into a kaon has a stricter upper bound on the respective branching ratio, 
namely 1.4 x 10~^, cf. Eq. (49). The solid turquoise line in Fig. 12 corresponds to this upper 
bound multiplied by a correction factor analogous to Eq. (46). We have taken into account that 
the BELLE collaboration required the kaon momentum to lie within a window of 1.6 GeV and 
2.5 GeV [74]. We can therefore exclude the MSSM parameter space above the turquoise line. 
Again, with a future super B factory we will be able to explore large regions of the parameter 
space in Fig. 12 [108]. 

To conclude this section, let us briefly comment on the validity of the mass insertion approxima- 
tion, Eq. (39). For the decay K~ — ^ T^'XiXiy the relevant mass insertion parameters are of the 
order 6i = O{10-^) and an expansion in 6^^ converges very fast. For the B decays, however, 
we have considered mass insertions of up to 5^^/^^ = 0.9, cf. Figs. 11 and 12. While it is clear 
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that the analysis should eventually be refined by taking into account higher order terms in the 
expansion Eq. (39), our approximation is sufficient to demonstrate that B decays into very light 
Xi can be tested with the next generation of B factories. 

V. LIGHT NEUTRALINO PRODUCTION AT HADRON COLLIDERS 

A massless or very light Xi that is purely bino is not fundamentally different in terms of collider 
signals than scenarios where it is (9(100 GeV) and mainly bino, such as SPSla (97 GeV, 98.5% 
bino). The main effect is that the phase space is bigger, in general increasing cross-sections, 
especially where the energy is not much greater than the creation threshold. There are many 
signals of interest at colliders, not least those involving decay chains, but the masslessness of 
the lightest neutralino does not fundamentally change those analyses. 

It is reasonable to ask whether the increase due to the greater phase space is in conffict with 
searches at colliders so far. The bino component does not couple to the Z boson at tree level and 
so constraints from the invisible width of the Z are only applicable to the very small higgsino 
components and also to the bino component at one-loop level. No bounds on the neutralino 
mass can thus be deduced from the Z width [16, 18, 20]. Direct searches at past and future 
e~^e~ colliders have been studied in detail in Refs. [20, 27, 28, 30, 32, 109]. While no bounds on 
the mass of the lightest neutralino can be obtained from past searches at LEP and B factories, 
measurements at a future e~^e~ linear collider might be able to discover light neutralinos through 
radiative production, e^e~^XiXi + 7- 

Direct hadroproduction of light neutralinos in association with a jet, 

pp/pp-^XiXi+ iet, (51) 

will lead to a spectacular monojet signature. We have calculated the tree-level cross-section for 
the pair production of massless Xi plus one jet, Eq. (51), for both the Tevatron {^/s = 1.96 
TeV) and the LHC {^/s = 14 TeV) using MadGraph [110]. We require a jet with transverse 
momentum of at least 80 GeV. The numerical results are collected in Table IV for the pseudo- 
SPS points described in Sect. Ill and are compared to the original SPS points with x5 masses 
of typically (9(100 GeV). Even though the Tevatron cross sections for massless x? are enhanced 
by a factor up to 2.5 with respect to the original SPS cross-sections with massive neutralinos, 
the expected number of events is only about 2 or less for the full 6 fb~^ of integrated luminosity 
of Run II so far. Given the large SM backgrounds, such as pp^Z^^uu) + jet, we conclude that 
current and future monojet searches at the Tevatron [111, 112] will not be sensitive to the direct 
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pair-production of light neutralinos. Similar conclusions hold for radiative production with an 
additional photon, pp^XiXi + 7- 

At the LHC, light neutralino pair production with jets is only enhanced from the massive 
case by about 20%. Still, with cross-sections of (9(100 fb), see Table IV, detection of these 
processes should be possible with sufficiently high luminosity and an excellent understanding of 
SM backgrounds [113, 114]. 



(pseudo-) 


Tevatron cross-section 


LHC cross-section 


SPS 


pseudo-SPS 


normal SPS 


pseudo-SPS 


normal SPS 


la 


2.96x10"^ fb 


1.68x10"^ fb 


3.26x10+2 fb 


2.73x10+2 fb 


2 


2.88x10"^ fb 


1.99x10-3 fb 


3.25x10+2 fb 


3.01x10+2 fb 


3 


2.66xl0~2 fb 


7.75x10-3 fb 


5.52x10+^ fb 


4.61x10+1 fb 


4 


4.48x10-2 fb 


1.84x10-2 fb 


8.59x10+1 fb 


7.14x10+1 fb 


5 


9.18x10^2 fb 


3.78x10-2 fb 


1.54x10+2 fb 


1.32x10+2 fb 



TABLE IV: Monojet cross-sections (T{pp/pp^XiXi + j^t) at the Tevatron and at the LHC. Shown are 
results for the pseudo-SPS points with massless Xi in comparison to the original SPS cross sections 
with massive neutralinos. 



VI. CONCLUSION 

Rare meson decays provide a sensitive test for the presence of new physics. In this paper we have 
studied the decay of pseudoscalar and vector mesons to light neutralinos in the MSSM. We have 
presented details of the calculations and explicit formulae for the two-body decays M— ^x^Xi^ 
where M = 7r°, rj, i]', p^, u, 0, J/ip, T. Furthermore, we have performed the first complete calcu- 
lation of the loop-induced decays K~^n~XiXi ^"^^ B''^K~XiXi- Considering various MSSM 
scenarios we find that the supersymmetric branching ratios are several orders of magnitude 
smaller then the SM processes with neutrinos instead of neutralinos in the final state. Con- 
sequently, no bounds on the neutralino mass can be inferred from rare meson decays in the 
MSSM with minimal fiavor violation. However, the branching ratios for the K~^tt~XiXi 
B^^K^XiXi decays may be significantly enhanced when one allows for non-minimal fiavor 
violation. We find new constraints on the MSSM parameter space for such scenarios and argue 
that future experiments in the kaon and B meson sector may be able probe light neutralinos 
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with masses up to approximately 1 GeV from rare decays. We have also considered searches for 
light neutralinos from monojet signatures at hadron colliders. While current and future monojet 
searches at the Tevatron will not be sensitive to the direct pair-production of light neutralinos, 
the detection of these processes should be possible at the LHC with sufficiently high luminosity 
and an excellent understanding of SM backgrounds. 
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APPENDIX A: FORM FACTORS FOR K' 7r-x?Xi AND 5" ^ K-/tt-XiXi 

For the decay K~ t^ XiXi^ the quark currents need to be replaced by their corresponding 
hadronic matrix elements involving a kaon K~ and a pion vr" [9] : 

{iT-{p^)\drPL/Rs\K-{pK)) = \ [umPK+p.r + f-mPK~p.r] , (ai) 

and 

/-W = [/o(t)-/+(t)]Mz:^, (A2) 

where t is defined via t = {px —PnY- The functions /+ and /q are the so-called Ki^ form factors. 
We describe them with the help of the linear approximation [9] 

/+/o(t) = /+/o(0) (^1 + A+/o^) , (A3) 

with A+ = 2.96 x IQ-^, Aq = 1.96 x lO'^ and /+(0) = /o(0) = 1.013, where we calculated /+(0) 
using the measured branching ratio for K~ — > 7T^e~Ue and isospin symmetry, i.e. 

{7r-{p^)\dYPL/Rs\K-{pK)) = V2{7t\p^)\uYPl/rs\K-{pk)) . (A4) 
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The decays 

B- ^n-/K-x',Xi, (A5) 

can be described in a similar way. One has to replace in Eq. (Al) px by and, if one is 
considering B^^K^x^Xi^ Pit by px- For the parametrization of the respective form factors 
f^^^ and fo^^ , we use [85, 86] 

r - ^ + 



= !i + 



K 



n/K 



fj^ = , (A6) 



where m^^ = 5.32 GeV and mf = 5.41 GeV are the masses of the B* and -B* vector mesons, 
respectively. The fit parameters are r[ = 0.744, = -0.486, {rrqf = 40.73 GeV^ rf = 0.162, 
rf = 0.173, rj = 0.258, rf = 0.330, (m^)^ = 33.81 GeV^ and (mf )2 = 37.46 GeV^. 
Using this parameterization for the non-MFV K~-^7t~XiXi decay, we obtain for the squared 
matrix elements summed over final state spins 



\Mll/rr\' = Cj^l/rr '' {^1+ + ^l- + ^--) ' (A7) 



with 



Ml+ = Afl[m'-o{m'^o + mi) + {h ■ A;2)(m^o - mi) 
+2m|o(A;i ■ p^) + 2m|o(A;2 ■ P-k) 
+2{ki ■Pn){k2 -Pjt)] , 



Ml^ = 8/+/_mio [m|o + ■ /cz) + {h ■ P.) + ■ p^)] 



Ml^ = 4/![m|o(A;i ■ k^) + m|o] , (A8) 

and the constants Cll = 1/1296 and Cjir = 1/81. The large difference between Cll and Crr 
originates from the different hypercharges of right- and left-handed squarks, with which they 
couple to the bino-like Xi- 
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